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Abstract The generation of harmonics by atoms in- 
teracting with two circularly polarized and frequency re- 
lated laser fields is addressed through ab initio numerical 
simulations. A detailed charaterization of a few specific 
harmonics is given. In particular, the two different cases 
where the total energy absorbed through photons is far- 
off or close to the energy gap between different atomic 
states are investigated. It is found that the conversion 
efficiency in the harmonic generation is strongly depen- 
dent on the inner atomic structure and in certain specific 
cases it can be significantly enhanced within a small fre- 
quency range. 



higher is i] the less is the number of harmonics in a cer- 
tain frequency range, and the circular polarization of the 
emitted harmonics. This last point becomes particular 
appealing for harmonics in the soft X-ray regime. Har- 
monic generation from sources other than atoms and cir- 
cular molecules, like linear molecules 0, nanotubes 
and plasmas |7j have been also investigated. The present 
paper is organized as follows: after a short summary of 
the theoretical frame within which our system lies, a 
full description of the model and of the numerical sim- 
ulations performed is presented. Finally, a conclusion is 
given. Atomic units (a.u.) are used throughout all the 
paper. 
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Introduction 

The interaction between an atom and two coplanar laser 
fields with circular polarization and commensurate fre- 
quencies is certainly an intriguing topic. In fact, it is 
worth to recall that when a single circularly polarized 
laser interacts with an atom no harmonic generation is 
possible. 

In the case where the two laser fields have opposite po- 
larization the emitted harmonics are of order 



n = k{ri + 1) ± 1, 



(1) 



where 77 is the frequency ratio and k G J\f^. Selection 
rules of this type are characteristic of all systems whose 
Hamiltonian is invariant under a dynamical symmetry 
operationjl], i.e., a transformation in both space and 
time. A further example of such systems is a circular 
molecule, e.g. benzene, interacting with a single circu- 
larly polarized laser field^jzlEMl- The derivation of such 
selection rules can be achieved by means of different ar- 
guments: angular momentum conservation, group the- 
ory and change of reference framed . The configuration 
we investigate here is of particular interest because of 
the filtering effect due to the selection rules i.e., the 



Theoretical frame 

In order to fully understand the nature of the hamon- 
ics we investigate, let us firstly summarize the deriva- 
tion of the selection rules ^ using angular momen- 
tum conservation. Let z be the propagation direction 
of the two lasers and ct^, the polarizations of the 
laser of frequency uj and rjuj, respectively. If a harmonic 
(7+ is emitted the sum of all the components along z of 
the angular momentum carried by the absorbed photons 
has to be -1-1. If the atom for example absorbs p pho- 
tons from the low-frequency laser, it must absorb p — 1 
from the other laser. Hence, the total absorbed energy is 
puj+ {p~ 1)770; — uj[p{l -t- 77) — 77] = uj[{p — 1)(1 + 1]) + 1]. 
Since p is arbitrary, with k = p + 1 and TV = 7; + 1 
we see that the harmonic of order {kN -I- 1) is emitted. 
With the same argument, starting from a harmonic of 
polarization ct^, we obtain kN — 1. The key point that 
should be stressed is the following: harmonic generation 
is possible if, and only if, the absorption of p photons 
from one laser is accompanied by the absorption of p ± 1 
photons from the other laser. 

Angular momentum considerations, as well as symmetry 
invariances, provide informations on the allowed orders, 
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but they cannot be used to derive any details concerning 
relative intensities and structures in the harmonic spec- 
tra like plateaus and cut-offs. Above all, stating that a 
certain order is allowed does not mean that the corre- 
sponding harmonic is really emitted. In fact, the previous 
arguments do not include case-specific features such as 
the atomic structure and the laser intensities. 
Recently, the characteristics and the field dependence of 
the harmonics emitted in our configuration have been 
addressed analytically through a calculation |S] based on 
the Lewenstein model Pl. With this analytical tool it is 
possible to derive a generalized cut-off for the harmonic 
spectra, an integral expression for the harmonic dipole 
strength, and to obtain again, independently, the selec- 
tion rules For what concerns the harmonic intensi- 
ties such a model treatment holds as long as the total 
absorbed energy by the electron is not "too close" to 
a transition[H|. In case the absorbed energy is far from 
any resonance different predictions obtained from the 
analytical model are confirmed through the numerical 
simulations. In particular it was found that if Ei and 
E2 are the electric fields of the low- and high-frequency 
laser, the intensity of the two harmonics with the same 
k are very different when Ei <^ E2, and they converge 
for El > £;2. 



Numerical simulations 

The numerical simulations were performed through an 
ab initio integration of the time-dependent Schrodinger 
equation on a two-dimensional (2D) grid. Choosing a 2D 
grid permits to run simulations quickly on every modern 
PC and is not expected to introduce qualitative modifi- 
cations to the phenomena we aim to investigate. In par- 
ticular, this holds for the specific configuration we are 
addressing here. In fact, both the electric field vectors of 
the lasers as well as the electron's trajectories lie in the 
plane that is represented by the 2D numerical grid. In the 
direction perpendicular to that plane the electron wave 
packet merely spreads. This spreading is not expected to 
strongly affect the main features of the harmonic gener- 
ation spectra such as the cut-off, for instance. Using a 
polar grid rather than a Cartesian one ensures that no 
artificial symmetry violation due to the numerics is in- 
troduced. 

In polar coordinates (p, </?), dipole approximation, and 
length gauge the time-dependent Schrodinger equation 
of our system reads 
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Fig. 1 Example of a typical harmonics spectrum obtained 
with two circularly polarized laser fields. The electric field 
intensities are E\ = E2 = 0.13. The frequency ratio is = 4 
and uj — 0.02 857r. No clear plateau structure is visible for 
these laser parameters. 



where the two laser pulses have a duration T — n/0, 
with ~ tj/84, and a sine-square shape. Vat{p) is a 
"soft-core" 2D potential given by 



Kt(p) = - 



(3) 



^{p,ip,t) 



(2) 



The parameters a and /3 can be tuned in order to ad- 
just the ionization energy. In our simulations we used 
a — 2.44 and /3 — 0.20. These values provide an ioniza- 
tion potential of Ip = 2.0, i.e., the one of real He"*". As 
we aim to address the role played by resonances in the 
conversion efficiency it is important to know the level 
scheme of the model potential we use. With the chosen 
parameters the lowest four excited states have energies 
f2i = 0.985, f22 = 1.375, = 1.548, and f2i = 1.592. 
These energies are measured with respect to the ground 
state. The laser frequencies have been chosen in accor- 
dance to this levels scheme. Concerning the frequency 
ratio rj all the results presented hereafter have been ob- 
tained taking 77 = 4. In Fig.n]an example of a two-color 
spectrum obtained with laser fields of equal intensities is 
shown. As expected only orders allowed by the selection 
rules are present. Substructures and broadening of cer- 
tain harmonics can be present when bound states other 
than the ground state come into playl*. 
The laser frequency ui has to be chosen carefully. In fact, 
on one hand we want the lower harmonics not to be af- 
fected by the atomic levels, on the other hand we want 
to approach resonances for the higher harmonics. The 
total absorbed energy which leads to harmonic genera- 
tion is LoN* with N* — k + ri{k±l). Obviously N* is also 
the order of the emitted harmonic. Therefore choosing 
oj — fii/N* provides the "virtual" resonance, i.e., the 
7V*-th harmonic is resonant with the transition between 
the ground state and the i-th excited state. This is a par- 
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Fig. 2 Energy scheme for the generation of harmonics no. 9, 
no. 11 and no. 14. The fundamental laser frequency is chosen 
so that the absorption process leading to the generation of 
harmonic no. 11 is resonant with the transition between the 
ground state and the field-free first excited state. 
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Fig. 3 Dipole strength for the harmonic no. 6 for different 
values of i?2. The laser frequency is ci; = 0.02 857r. The behav- 
ior is purely perturbative. 

ticular kind of resonance, in fact none of the two absorp- 
tion processes is directly resonant, but their combination 
is. The ionization rate is nearly not affected by the close- 
ness of such kind of resonances and remains always very 
low. Choosing N* = 11 and i = 1 the expected resonant 
frequency is w = 0.02857r. A scheme of such a configura- 
tion is presented in Fig.|21 While the harmonic no. 11 is 
exactly resonant with the first excited state, harmonics 
no. 6 and no. 9 are not affected by any energy level. Once 
the frequency value was fixed, a series of simulations for 
different laser fields intensities have been performed. As 
we deal with two different absorption processes, in or- 
der to obtain efficient harmonic generation it is neces- 
sary that the two separate processes are "likely" , i.e., the 
probabilties of absorbing the required numbers of pho- 
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Fig. 4 Dipole strength for the harmonic no. 11 for low (a) 
and high (b) values of E^ - The laser frequency is w = 0.02 857r. 
A perturbative behavior is observed only in the case of rela- 
tively weak fields. 



tons from each of the two laser must be of the same order 
of magnitude. 

In Fig.|31 the behavior of the harmonic no. 6 versus 
the electric field E\ for different values of Ei is plot- 
ted. As expected the intensity of the harmonic increases 
smoothly with increasing laser intensity. A very differ- 
ent behavior is observed for the harmonic no. 11 (Fig.01). 
For relatively low intensities of E\ the behavior is simi- 
lar to that of the harmonic no. 6, but when the electric 
field E\ becomes stronger a much more complex behav- 
ior appears. This is clearly due to the resonance with the 
first excited state. Actually, the frequency lo = i7i/ll is 
close to resonance, but not exactly resonant. In fact, one 
should take into account the level shifts due to the dy- 
namical Stark effect: all levels are expected to be moved 
upward. Giving an analytical estimation on how much 
the lower lying excited states are moved up is a very dif- 
ficult task already when only a single, linearly polarized 
laser field is taken into account, let alone in our two-color 
configuration. In particular it is difficult to distinguish 
between the different contributions of the two laser fields. 
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Fig. 5 Dipole strength of harmonic no. 11 vs i5i for differ- 
ent laser frequencies. E2 = 0.13 is held constant. A small 
increase (or decrease) of the laser frequency affects strongly 
the efficiency conversion of the harmonic which approaches 
the resonance. 
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Fig. 6 Dipole strength of harmonics no. 9 and no. 11 vs laser 
frequency. A strong enhancement of harmonic no. 11 due to 
the resonance is observed. 

Conclusions 



A possible method for studying how a full resonance af- 
fects the harmonic generation is to perform a series of 
numerical simulations where all parameters but the laser 
frequency w is kept constant. Increasing the energy of 
the photons compensates the shift due to the Stark ef- 
fect and allows to reach the shifted excited state again 
with 11 photons. As the amount of the shift is unkwown 
different values of w have been used. The results of these 
simulations are presented in Fig.[S| It appears that the 
harmonic intensity strongly depends on the exact value 
of the laser frequency, i.e., a small increase or decrease 
in the frequency value can change the harmonic inten- 
sity several orders of magnitude. From Fig.El one real- 
izes that the full resonance is achieved for ll> — O.OSOtt 
and E\ = 0.07, leading to an energy shift which is about 
0.052 a.u. All these values are obtained keeping constant 
E2 = 0.13. If the value of E2 is different the previous 
values for uj and Ei do not hold anymore although the 
physics remains qualitatively the same. In order to bet- 
ter estimate the dependence of the conversion efficiency 
with respect to the exact location of the excited state 
a different kind of study has been carried out. Taking 
the laser intensities which in Fig.|31 give the maximum 
harmonic intensity, namely Ei = 0.07 and E2 = 0.13, a 
series of simulations for different laser frequencies within 
a small frequency range has been performed. The results 
of such study are presented in Fig.El While the inten- 
sity of harmonic no. 9 varies slowly, the intensity of har- 
monic no. 11 exhibits a strong enhancement due to the 
resonance. The ionization level remains negligible for all 
frequencies values presented in Fig.|51 



In this work an extensive numerical study of harmonic 
generation in the two-color coplanar configuration has 
been presented. In particular the cases of far-off and near 
resonances absorption have been addressed. It has been 
shown that the atomic levels can be used as an important 
tool in order to enhance significantly the intensity of 
a particular harmonic without increasing the ionization 
rate. 
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